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Abstract 

We prove a rigorous upper bound for the effective conductivity of an isotropic composite made 
of several isotropic components in any dimension. This upper bound coincides with the Hashin 
Shtrikman bound when the volume ratio of all phases but any two vanish. 

1 Introduction 

The effective conductivity of a composite material made of two isotropic phases with given conduc- 
tivities and volume ratios is bounded below and above by the Hashin-Shtrikman bounds [3]. These 
bounds are exact. For any value of the parameters there is some micro structure that realizes the 
bounds. The Hashin-Shtrikman bounds can be extended to composites of more than two phases, 
but they are not exact for some values of the parameters. As it was pointed out by Milton [9], the 
Hashin-Shtrikman upper bound for three phase composites for which the most composite phase has 
volume ratio zero, does not coincide with the bound for the two other phases only. This is highly 
counterintuitive, and strongly suggest that there should be better bounds for some values of the 
parameters. 

In two dimensions, the Hashin-Shtrikman bounds for multiphase composites were refined by Nesi 
[10] and Astala and Nesi [JJ . Exact bounds for the two dimensional case were proposed by Cherkaev 
in [2]. In more than two dimensions, there is no bound that refines the Hashin-Shtrikman ones. The 
purpose of this paper is to progress in that direction. 

In this paper we prove a new upper bound for the effective conductivity of an isotropic composite 
with an arbitrary number of isotropic phases in arbitrary dimension that refines the Hashin-Shtrikman 
bound. The bound is not exact except when Hashin-Shtrikman bounds are exact. In fact the bound 
seems to be very rough, but it does satisfy the right asymptotics. If we let the volume ratio of all 
phases but two go to zero, the bound converges to the corresponding two-phase Hashin-Shtrikman, 
for any choice of the two phases to keep. This provides the first rigorous proof, in dimension larger 
than two, of the intuitively obvious fact that a phase with a negligible volume ratio has a negligible 
effect in the possible effective conductivities. 

We consider periodic composite materials made from isotropic phases. The microscopic structure 
is given by a periodic function a. Typically, for a composite with K phases, the function a has the 
following structure. 

<j\ in A\ 

a(x)=( (1.1) 
k a K in A K 

The sets A%, ■ ■ ■ , Ak are a partition of the unit cube Q of K™ and <r(x) is extended periodically 
outside Q. We denote by fii the measure of each set Ai. Obviously, since {J i Ai = Q,/j,i + -- -+HK = 1- 



1 



The effective conductivity associated with a is computed via the celi problem: 



(Av,v) = min / a(x)\v + Vitr dx (1.2) 
ueHl(Q)J Q 

where Hp{Q) stands for the H 1 functions in the unit cube Q with periodic boundary conditions. 

In this article we will concentrate in finding bounds for the trace of A that depend on the values 
of [ii and <Ji but not on the particular structure A,. In case the composite has cubic symmetry, 
then the effective conductivity is isotropic. In that case, our bounds are just bounds of the effective 
conductivity. 

We define 

1 . 1 



o=-XxA=~ \ (Ae t ,e t ). 
n n . ^— ' 



Our purpose is to find upper and lower bounds for a. 

Clearly, the multiphase composite problem is a particular case of the problem of finding bounds 
for a when all that is known of the periodic coefficient a is its distribution function 

F(t) = \{x : a{x) > t}\. 

Thus, in this work we will consider the problem of finding an upper bound of a that depends on F 
only. It is important to remember that inf a, sup a and an integral of the form 

G(a(x)) dx, 

for any function G, depend on the distribution of a only. 

The following theorem provides an upper bound that refines the Hashin-Shtrikman bounds. 

Theorem 1.1. For any value of the parameter S > 0, the following upper bound holds: 

a <H + E 

where 



H = -(n - 1)5 + (J (a(x) + (n - 1)5)' 

C(osca) 2 (f (a(x) + (n-l)S)- 1 )~ A r 

E = c f liw T7 4 / - ^gF(t)f dt 

(inf a + (n— l)A)^(supcr + (n — l)o) Js 

where C is a constant depending only on dimension. 

We will not provide an explicit expression for the constant C in Theorem ll.il It is a constant that 
depends on dimension only. It is related to the BMO estimates for the Laplace equation in the unit 
cube with a bounded right hand side. The best constant for these estimates is difficult to compute. 

The proof of Theorem ll.il uses the idea (originally from [11]) that an upper bound can be obtained 
by reducing the set of possible vector fields in the cell problem to those that are the gradient of a scalar 
potential function. The term E in the upper bound of Theorem 11.11 is obtained applying harmonic 
analysis estimates to this potential function. With this approach, upper bounds are more natural to 
obtain than lower bounds (unlike other approaches for bounds like the translation method). In this 
article we leave the question of a similar lower bound open. 
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Note that for any S < sup a, 



(n - l)S < (^J (a(x) + (n - < (n - 1)5" + supo-. 



Therefore we can simplify the expression of E for the bound in theorem 11.11 to 

C(ose g ) 2 r Fm -i ogFit) f dt (1.3) 



(infa+(n-l)5) s 



s 



Let us see what bounds Theorem 11.11 provides in the case of a three phase composite. In this case 
the function a has the form 

{o\ in A\ 
02 in A 2 
a 3 in A 3 

where \Ai\ = /xi, |A 2 = \ii and A 3 = /Lt 3 . 

If we choose S = 03, then the term E in Theorem 11.11 equals to zero, and we recover the classical 
Hashin-Shtrikman bound for three phase composites. 

If we choose S = 02 , then we obtain the following bound using the upper bound (|1.3[) , 

. , , , -n i-i I , C(cr 3 - <7i) 2 (a 3 - CT 2 )/i 3 (l - logM3) 2 

a < -(n- 1)0-2 + > /i. d.i n-lH + ] — — 7 ■ ( L4 

I ' ' (cti + (n - l)cr 2 )^ 



. i=l,2,3 



This bound refines the Hashin-Shtrikman bound when the volume ratio /z 3 is small. Indeed, as 
/i 3 — > 0, (|1.4j) converges to the two phase Hashin-Shtrikman bound corresponding to the first and 
second phase only. 

The organization of the paper is as follows. In section[2]we prove the trivial upper bound that a is 
bounded above by the average of a. In section[3] we provide a proof of the classical Hashin-Shtrikman 
bounds, and in section|3]we prove the new bounds of Theorem ll.il Even though the bounds obtained 
in sections [5] and [3] are well known, we use these two sections to introduce the ideas that will be used 
to obtain the new bounds in section 01 These ideas originated in [TT]. The inequalities obtained in 
section [3] are the basis of the proof of Theorem 11.11 in section [4] 

2 The trivial upper bound 

Let us analyze the formula for a . 

a = — (Aet, ei) 

= — min / <j(x)\ei + Vud 2 dx 

= - min / a(x)\I + DU\ 2 dx 
nu&Hi{Q,9^)J Q 

were Hp(Q,M. n ) stands for the space of periodic function is Q with values in K n and finite norm in 
H . Indeed U = (ui, . . . , u n ). Therefore, the quantity that we have to estimate is 

a=- min f a(x)\I + DU\ 2 dx (2.1) 
n ueHi(Q.R")J Q 
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One way to obtain an upper bound for (|2.1[) is by using a sample vector field U. The trivial upper 
bound is what we obtain by setting [7 = 0. 



a < / <j(x) Ax 
Jo 



In the case of a composite with K phases, this trivial upper bound becomes X)j=i k A 4 * *- 
This trivial upper bound is not achievable as soon as there is more than one different phase. 

3 The Hashin-Shtrikman upper bound 

We will start by showing a simple way to obtain the standard Hashin-Shtrikman bound for a . The 
idea of this proof was first given in . This idea was extended in [B] where it was used to obtain the 
full Hashin-Shtrikman bounds for multiphase composites and study their attainability conditions. 

The main idea is to obtain an upper bound for (|2.I[) by reducing the set where we look for a 
minimum to only those vector fields U that are the gradient of a potential p. Therefore 

a<- min / a(x)\I + D 2 p\ 2 dx. (3.1) 
n P eH}(Q)J Q 

The following elementary relations will be used. Note that the first identity plays a key role in 
the usual proof of the Hashin-Shtrikman bounds by the method of compensated compactness ([7], [5] 
and [T3] ■ See also [5]). 

• For any periodic function p in H 2 we have 

\D 2 p\ 2 dx= { \l\p\ 2 dx (3.2) 



Q JQ 

• For every matrix M , the following inequality holds 

\M\ 2 > -\trM\ 2 . (3.3) 
n 

In particular, for M = D 2 p, \D 2 p\ 2 > ^\Ap\ 2 . The equality holds only if D 2 p is a scalar matrix. 

• For any periodic function p 

[ Ap = (3.4) 
Jq 

Moreover, any function with zero average in the cube is the Laplacian of some potential p. 
Let S be an arbitrary constant. We apply these relations to (|3.1[) and obtain 

a=— min / a(x)\I + D 2 p\ 2 dx 
n P eH}(Q)J Q 



>Q 

min h(Ap) + I 2 (D 2 p) 
pem(Q) 



(3.5) 



where 



h(Ap) = - [ a{x)n + 2a{x)Ap + S\Ap\ 2 + -(cr(x) ~ S)\Ap\ 2 dx, (3.6) 
' n J Q n 

h{D 2 p) = - [ (a(x) - S) (\D 2 p\ 2 - -(Ap) 2 ) dx. (3.7) 
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The advantage of these estimate is that (|3.6[) depends only on Ap. The value of Ap(x) can be 
chosen arbitrarily for every point x with the only condition that it has average zero on Q. 

In order to find an upper bound for er, we will choose the p that minimizes ii(Ap), and then we 
will estimate the remainder I2 {D 2 p) using elliptic estimates for the obtained function p. 

Since Ap has average zero on Q, we can rewrite I\ as 



I 1 {Ap) = - [ <T(x)n + 2(a(x) + (n-l)S)Ap + S\Ap\ 2 + -(a(x)-S)\Ap\ 2 dx, 

IQ " 

Ap x 2 



n 



= J a(x) + (a(x) + {n-l)S)Ul+^j -lj dx, 

= -(n - 1)5 + J (a(x) + (n - 1)5) + dx 

Now we will find the value of Ap that makes the value of I\ smallest. We compute the first 
variation of 1\ : 

Dh-h = J (<t(x) + (n - 1)5) ^1 + — J h(x) dx, 

which has to be zero for any function h with average zero on Q. This means that there is a constant 
L such that 

(a{x) + (n-l)S) (i+^f) = L - 
We can compute the exact value of L using again that Jq Ap = 0. 



IQ 

Now, for this optimal choice of Ap, we have 



= (| (a(x) + (n-l)S)-^ 1 



(3.8) 



h(Ap) = -(n-l)S + J L^l + ^j dx 



(3.9) 

= -(n - 1)S + L = -(n - 1)5 + l^J (a{x) + (n - 1)5) -1 ^ (3.10) 

Let us now concentrate on the value of I2 (D 2 p) for this same choice of function p that minimizes 
h(Ap). 

We first observe that from (|3 . 3[) . (|£> 2 p| 2 — i(Ap) 2 ) > 0. Therefore, we can disregard the points 
x where a(x) < S: 



h{D 2 p) < - I (a(x) - S)+ (\D 2 p\ 2 - -(Ap) 2 ) dx. 
n Jo V n ) 



(3.11) 



In particular if we choose S = sup a (= <7k in the case of a multi-phase composite) we make 
I 2 (D 2 p) — 0, we obtain that a < Ji(Ap) and we recover the Hashin-Shtrikman bound 



a < —{n — 1) super + (^J (cr(x) + (n — 1) super) ^ 



(3.12) 



The attainability condition is that the minimal vector field U in (|2.1[) is the gradient of a scalar 
potential p and also Z? 2 p is a constant scalar matrix in all phases except the most conductive one. The 
range of parameters for which one can guarantee the attainability of the bounds (|3.12[) is analyzed in 
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4 New bounds 



In the optimal configuration for the Hashin-Shtrikman bounds, the most conductive phase plays 
the role of the matrix and all other phases are inclusions. The Hashin-Shtrikman bound cannot be 
optimal for multiphase composites in all range of parameters. Indeed, if we apply (|3.12[) to a three 
phase composite and make the volume ratio of the most conductive phase ^3 approach zero, the 
bound (|3.12l) still depends on 03 in the limit, even though the third phase is not present. When the 
most conductive phase has a very small volume ratio, it is more convenient to choose a different value 
of the arbitrary parameter S and estimate (13.71) . 

In this section we prove theorem ll.il The proof refines the estimate from section [3] In particular 
it is based on the estimates (I3.5[) . (|3.6p and p. lip . In order to find an upper bound for Ii{D 2 p) we 
will need the following technical Lemma about BMO functions. 



Lemma 4.1. Let f : Q 



be a BMO function in Q with average zero. Then for any subset A C Q, 



I |/| a da:<C| 

J A 



\BMO 



(l-log|A|) 2 L4| 



(4.1) 



where \A\ stands for the measure of A and C depends only on dimension. 



Note that (|4.1j) also holds if / is vector or matrix valued, since the inequality can be applied 
component- wise. We defer the proof of this Lemma to the appendix. We now prove the main 
theorem. 

Proof of Theorem \l.l[ As we mentioned above, we will find an upper bound for a using (|3 . 5[) . Our 
upper bound corresponds to one test function p, which may or may not be optimal for (|3.5[) . In fact, 
we will choose the function p which is optimal only for the value of I\(Ap) and then we will estimate 
the value of l-iij^p) from above. 
Note that 

\D 2 p\ 2 - -\Ap\ 2 \ = \D 2 p - i(tr D 2 p)I\ 2 . 



Thus, l2{D 2 p) is the integral of the squared norm of the traceless part of D 2 p. 
Recall from section [31 that the potential p that minimizes I\(Ap) solves 

Ap = nL{a(x) + (n - 1)5) _1 - n =: 9{x) (4.2) 

For this p, it is simple to estimate the L 2 norm of D 2 p— -Apl using the Fourier transform. Indeed 



therefore 



D 2 p- - Apl 
n 



1 



(0 = - 



lei 2 



D z p- - Apl 



< 



L 2 



n- 1 



n 



1/2 



lam) 



\L 2 



An L 2 estimate for the integrand only provides a rough estimate for (I3.1ip that would only take 
into account the maximum value of (o~(x) — S). On the other hand, 9 is a bounded function. We can 
use Calderon-Zygmund theory to obtain that D 2 p is a BMO function whose norm depends on the 
oscillation of 9 only (see |12)b 



D 2 p Apl 

n 



BMO(Q) 



< C osc 9 
Q 



(4.3) 
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where C is a constant depending only on dimension. 
We compute the oscillation of 8 from (|4.2[) : 



osc 8 = nL 
Q 



1 



inf a + (n — l)S super + (n — l)S 

n 



(4.4) 



(inf a+(n-l)S) (sup a + (n - 1)5) 



LOSCCT 



Recall that L was given in (|3.8[) . 

Combining (|4.4|) with (|4.3|) and (|4.1|) we can obtain an upper bound for l2(D 2 u) that depends on 
the distribution function F of a only. From Fubini's theorem: 



h{D 2 p) < - f (a(x) - S)+ (\D 2 p\ 2 - -(Apf) dx, 



1 

< - 

n 

1 

< 



■ 



a(x)>t 

1 



D 2 p\ 2 - -(Ap) 2 ) dxdt 
n 



D z p Apl 

n 



BMO JS 

CL 2 ( osc a) 2 



F{t){\ ~ \og F {t)f dt 

" + 00 



(infer + (n - l)S') 2 (sup <r + (n - 1)5) 2 



F(t)(l-logF(t)) 2 dt 



Adding this bound for li(F) 2 p) with the bound (|3 . 10[) for Ii(Ap) we finish the proof of Theorem 
LU □ 



Remark 4.2. The constant C in (|4.3|) must be smaller for the traceless part of the Hessian: 

\\D 2 p — ^Apl\\ BMQ ^ than for the whole Hessian ||-D 2 p|| BMO (q)- It is hard to know the best 

constants though. 

Remark 4.3. The obtained bounds for l2(D 2 p) seem to be very crude. In (|4.3p we estimate the 
BMO norm of D 2 p with respect to the L°° bound of Ap. An estimate in L°° for D 2 p is known 
to be false for the Laplace equation. On the other hand, as it is pointed out in [2], if there exists a 
composite with piecewise differentiable interfaces where the bounds for (12 . 1[) are realized, \DU\ should 
be bounded. 



5 Appendix: The proof of Lemma 14.1 

In order to prove Lemma 14.11 we will use the following classical result by F. John and L. Niremberg 
about BMO functions (for a proof see [4]). 

Lemma 5.1. Let f be a BMO function onQ = [0, 1]™ with average zero. Then the following inequality 
holds for the distribution of f : 

\{\f \ > o-}\ < Be W^bmo 
where b and B are constants depending only on the dimension n. 
Proof of Lemma \4- 1\ Let T > be a parameter which will be chosen later. 
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We use Lemma I5TT1 to estimate the integral. 

/ |/| 2 dx= [°° 2t\{\f\ > t} n A\ dt 

J A JO 

< J 2tmin (Be imLawo , \A\j dt 

For any T G (0, +00), 

/* Z* 00 -bt 

/ |/| 2 da; < / 2t\A\dt+ / 2tBe^^BMo dt 

J A JO JT 

= T 2 \A\ + 2b(t + Mlmm j l^kMQ e TCi^ 

Choosing now T = - lml ^ MO log \A\, 

<C\\f\\ BMO 

(l-logL4|) 2 |A| 

for a constant C depending only on dimension. □ 
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